Thermodynamics of the superconducting state in Calcium at 200 GPa 
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The thermodynamic parameters of the superconducting state in Calcium under the pressure at 
200 GPa were calculated. The Coulomb pseudopotential values (fi*) from 0.1 to 0.3 were taken into 
consideration. It has been shown, that the specific heat's jump at the critical temperature and the 
thermodynamic critical field near zero Kelvin strongly decrease with fi* . The dimensionless ratios 
n = AC (T c ) /C N (T c ) and r 2 = T C C N (Tc) /H c (0) significantly differ from the predictions based 
on the BCS model. In particular, n decreases from 2.64 to 1.97 with the Coulomb pseudopotential; 
(3JT)' whereas r 2 increases from 0.140 to 0.157. The numerical results have been supplemented by the 

^ , analytical approach. 
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PACS numbers: 74.20. Fg, 74.25.Bt, 74.62. Fj 



I. INTRODUCTION 
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By using the advanced technique, it is experimentally possible to explore the properties of the superconducting 
Q_i' state under the high pressure (p). In particular, the above researches enable: (i) the test of the theories for the 
superconducting state, (ii) as well as to improve the properties of the superconductors, and (iii) to create the new 
superconductors. At present, the 52 elemental superconductors are known, however 22 of them superconduct if the 
pressure is applied [l[. The most interesting elements are: Lithium, where the critical temperature (Tc) rises rapidly 
to ~ 14 K at 30.2 GPa Yttrium with the maximum value of the critical temperature of 20 K at 115 GPa and 
Calcium which has the highest observed value of Tc (25 K at 161 GPa) [1], @, @. 
. The thermodynamic properties of the superconductors under the high pressure can be analyzed in the framework 
of the Eliashberg approach 0, H) @- In this formalism the complicated form of the electron-phonon interaction 
is modeled by Eliashberg function (a 2 F(fl)). We notice that the first-principle calculations of a 2 F (f2) require the 
i knowledge of the electronic wave functions, the phonon spectrum, and the electron-phonon matrix elements between 

two single-electron Bloch states. Experimentally, the form of the Eliashberg function can be directly obtained from 
the second derivative of I-V curve for the tunnel junction (d 2 V/dI 2 ~ a 2 F (mV)) [l(J. 

From the physical point of view, the Eliashberg approach represents particularly important method of the analysis, 
since it enables the calculation of the thermodynamic parameters on the quantitative level. In particular, the exact 
form of the free energy difference between the superconducting and normal state should be calculated on the basis 
of the so-called Eliashberg equations [ll[ . In the considered case, the input parameters are the Eliashberg function 
and the Coulomb pseudopotential (fi*), where fi* models the Coulomb repulsion between electrons. We notice that 
its value is selected in such way that Tc determined on the basis of the Eliashberg equations equals the experimental 
value of the critical temperature. 

In the presented paper, we have analysed the thermodynamic properties of Calcium under the pressure at 200 GPa 
• • . by using the Eliashberg approach. In particular, the following parameters were taken into account: the specific heat 
. m the superconducting state (C s ), the specific heat in the normal state (C N ), and the thermodynamic critical field 
(H c ). Additionally, the dimensionless ratios n = (C s (T c ) - C N (T c )) /C N (T c ) and r 2 = T C C N (T c ) /H 2 (0) have 
5_i ' been determined [12 ]. 

For Calcium, the dependence of the critical temperature on the pressure has been obtained experimentally by Okada 
et al. in 1996 Q and then by Yabuuchi et al. in 2006 Q. The results prove that the critical temperature grows with 
the pressure from the value of about 1 K (p = 50 GPa) to 25 K (p = 161 GPa). On the basis of the Yabuuchi's results 
it is easy to show that the values of the Coulomb pseudopotential can be large. For example, ]/ i *] p= i20GPa = 0.215 and 



[M*] D= i60GPa = 0.241 [13| . [I4j . In the case of the pressure 200 GPa, the experimental value of Tc is still unknown, 



therefore the wide range of the Coulomb pseudopotential's values have been considered in the paper; ^* g (0.1,0.3). 
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II. THE ELIASHBERG EQUATIONS 



The Eliashberg functions for Calcium were determined by Yin et al. @ . At the pressure 200 GPa, the Prima struc- 
ture is clearly favored, and the linear-response calculations indicate that it is also dynamically stable. Additionally, 
the Yin's results have shown that the strong electron-phonon coupling persists and Tq can be high (~ 30 K). In the 
paper, the thermodynamic parameters for Calcium have been calculated by using the Eliashberg equations on the 
imaginary axis, the Yin's Eliashberg function has been taken into account. 

The Eliashberg set has the following form: 
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P m =-M V^nZl + 0: 

The wave function renormalization factor is denoted by Z n = Z {iuj n ), and the order parameter function by <j> n = 
(j>(iu> n ). The n-th Matsubara frequency has the form: u n = (n/0)(2n— 1), where = (ksT) and fcs is the 
Boltzmann constant. We notice that the value of the order parameter is given by the ratio: A„ = <fi n /Z n . In Eqs. (Q]) 
and ((2]) the symbol A (z) denotes the electron-phonon pairing kernel: 

A(z) = 2 jT nm " dn-^-^a 2 F(n), (3) 

where the maximum phonon frequency il max is equal to 78.1 meV. Finally, 9 denotes the Heaviside unit function and 
lo c is the cut-off frequency (uj c — 3il max ). 

From the mathematical point of view the exact solution of the Eliashberg equations represents a complicated 
problem. We notice that formally the Eliashberg set contains the infinite number of the non-linear algebraic equations; 
in addition every equation has the integral kernel, which is dependent on the form of the Eliashberg function. It is 
possible to prove, that if we limit the number of the Matsubara frequencies, the solutions of the Eliashberg equations 
lose the convergence only in the area of the very low temperatures. In the paper we assume M = 1100. In this case 
the functions Z n and 4> n are stable for T > 3.48 K. The Eliashberg equations have been solved by using iterative 
method [151 ] . 



III. THE NUMERICAL RESULTS 



In Fig. [T]the order parameter as a function of the number m has been presented. We have considered the selected 
values of the Coulomb pseudopotential and temperatures. It is easy to see that A m decreases with m number's growth. 
In particular, the strong fall of the order parameter appears for low values of m; for higher values the function A m 
saturates. Very similar dependence on the number m possesses the wave function renormalization factor (see Fig. 
However, the function Z m saturates considerably slower than A m . 

The growth of the Coulomb pseudopotential's value differently influences on the order parameter and the wave 
function renormalization factor. In the first case the values of the order parameter strongly decrease; whereas the 
function Z m is practically not changing. From the physical point of view the above facts mean, that together with 
the increasing of fi* decreases only the value of the critical temperature; the electron effective mass remains fixed 
(m* - Z m=1 ). 

The solutions of the Eliashberg equations also very unlikely evolve with the temperature. Most clearly this fact is 
possible to observe in Fig. [3l where the dependences of A m= i and Z m —\ on the temperature have been plotted. In 
particular, the obtained results show that the temperature dependence of the order parameter can be modeled by the 

function: A m= i (T) = A m= i (0) yl~ (i^) ' wnere tne parameters A m= i (0), Tc, and are collected in Tab. Q] On 

the other hand, the wave function renormalization factor is slightly depended on the temperature. 

In order to obtain the specific heats and the thermodynamic critical field, we have to calculate the free energy 
difference between the superconducting and normal state (AF) [ll[ : 



AF 
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FIG. 1: The dependence of the order parameter on m for selected values of the Coulomb pseudopotential and temperatures. 




FIG. 2: The dependence of the wave function renormalization factor on m for selected values of the Coulomb pseudopotential 
and temperatures. 



The functions and Z% represent the wave function renormalization factors for the superconducting (S) and normal 
(N) state respectively. 

The specific heat difference between the superconducting and normal state (AC = C s — C N ) can be obtained by 



TABLE I: The values of A m=1 (0), T c , and f3 parameters. 
A m =i (0) mcV T C K 
0.1 6.48 36.15 3.47 

0.2 4.63 26.94 3.50 

0.3 3.45 20.79 3.60 
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FIG. 3: The dependence of the order parameter and the wave function renormalization factor on the temperature for selected 
values of the Coulomb pseudopotential. In both cases we assume m — 1. 



using the expression: 



AC 



Id 2 [AF/p(0)] 
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On the other hand, the specific heat in normal state is given as 
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where 7 = |7r 2 (1 + A). In Fig. |4]the specific heat for the superconducting and normal state as a function of the 
temperature has been shown. It is easy to see that the specific heat's jump at the critical temperature decreases with 



AC(T C ) 
P(0) 



/ 



the growth of the Coulomb pseudopotential. In particular 

^i*=0.3 

The thermodynamic critical field has been calculated by using the formula 



AC(Tc) 
P(0) 



0.43. 



VpW) 



V-Stt [AF/p(0)]. 



(7) 



In Fig. [5] the dependence of Hq/ \J p (0) on the temperature has been presented. We can see, that the value of the 
thermodynamic critical field near the temperature of zero Kelvin (Hq (0) ~ Hq (To)) also strongly decreases with /x*; 



/ 



=0.3 



0.55. 



Ai*=0.1 



The dimensionless ratios r\ and on the basis of the calculated thermodynamic functions have been determined. 



We notice that in the framework of the BCS model these parameters have the universal values: [ri] 



BCS 



1.43 and 



[?'2] BCS = 0.168 [12|. For Calcium the theoretical data have been collected in Fig. [6l We see that even for large values 
of [i* the ratios significantly diverge from the values predicted by the BCS model. Below, we have given the formulas, 
which enable the simple calculations of r\ and r^. In particular: 



= 1-42.5 
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FIG. 4: The dependence of the specific heat in the superconducting and normal state on the temperature for selected values of 
the Coulomb pseudopotential. The vertical line indicates the position of specific heat's jump that occurs at Tc- 
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FIG. 5: The thermodynamic critical field as a function of the temperature for selected values of the Coulomb pseudopotential. 



and 



1 + 2.7 



BCS 



II Win 
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(9) 



We notice that the numerical coefficients in Eqs. © and © by the least-squares analysis of 60 exact n and 
ri values have been chosen. Finally, the critical temperature by using the modified Allen-Dynes formula should be 
calculated fill: 



ksTc = /i/2 7^exp 
1.45 



-1.03(1 + A) 



A-^*(1 + 0.06A) 



(10) 
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FIG. 6: The ratios r\ and V2 as a functions of the Coulomb pseudopotential. The black circles show the exact results obtained 
on the basis of the Eliashberg equations. The solid lines represent the calculation of the n and ri parameters using our 
analytical scheme (Eq. ((8} and Eq. (JSJ). The dashed lines represent the results obtained in the framework of the Marsiglio [ij 
and Carbotte [t| formulas respectively. The dotted lines denote the BCS values. 



where the strong-coupling correction function (/i) and the shape correction function (/ 2 ) are given by: 



and 



The functions Ai and A2 have the form: 



and 



(11) 



Win 



1 X 



A 2 + A2 



Ax = 2.6 (1 + 1.8//*), 



A 2 = 0.092 (1 - 150/**) 



Win 



(12) 



(13) 



(14) 



The parameter L02 is the second moment of the normalized weight function, u>i a denotes the logarithmic phonon 
frequency and A is called the electron-phonon coupling constant. In the case of Calcium under the pressure at 200 
GPa the following results have been obtained: y/Zj-J = 35.92 meV, u)\ n = 29.98 meV and A = 1.228. 



IV. SUMMARY 



The thermodynamic parameters of the superconducting state in Calcium under the pressure at 200 GPa have been 
analyzed in the paper. The numerical calculations in the framework of the Eliashberg approach have been made. 
On the basis of the exact Eliashberg solutions the specific heats and the thermodynamic critical field have been 
determined. 
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For the wide range of the Coulomb pseudopotential values (fi* € (0.1,03)), it has been shown that the specific 
heat's jump at the critical temperature and the thermodynamic critical field near zero Kelvin strongly decrease with 
/'•• 

The dimensionless ratios t\ and ri significantly differ from the predictions based on the BCS model even for high 
values of the Coulomb pseudopotential. In particular, t\ g (2.64, 1.97) and r 2 € (0.140,0.157). 
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